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■^j- , Abstract 

Recently, B. Gerganov, A. LeClair and M. Moriconi [Phys. Rev. Lett. 86 (2001) 4753] have proposed an "exact" 
(all-orders) /3-function for 2-dimensional conformal field theories with Kac -Moody current-algebra symmetry at any 
level k, based on a Lie group Q, which are perturbed by a current-current interaction. This theory is also known as 
the Non-Abelian Thirring model. We check this conjecture with an explicit calculation of the /3-function to 4-loop 
order, for the classical groups Q = S\J(N), SO(iV) and SP(iV) at level k = 0. We find a contribution at 4-loop order, 
proportional to a higher-order group-theoretical invariant, which is incompatible with the proposed /5-function in all 
possible regularization schemes. 
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1 Introduction 



Perturbations of conformal field theories (CFT) in two dimensions (2D) have been a very active topic of study for 
a long time. The focus of this paper are 2D conformal field theories possessing Kac-Moody current algebra (or: 
"Affine Lie algebra") symmetry [1] associated with a Lie group Q, which are perturbed by a bilinear in the Noether- 
current (i.e. by a "left-right current-current bilinear interaction"). Non-abelian Thirring models [2,3] and Gross- 
Neveu models [4] are much studied examples. Typically, such perturbations are (marginally) relevant and generate 
a mass scale; in these cases the long-distance (infrared) behavior of these theories is that of a massive field theory. 
However, generalizations of these theories may exhibit [5] massless long-distance fixed points in a certain zero 
species ("replica") limit, or, which can be seen to be equivalent [6], when the symmetry group Q is replaced by 
a certain supergroup. Such theories are of great interest in condensed matter physics, since the mentioned "zero- 
species", or equivalently, supersymmetry generalizations describe (de-)localization transitions known to occur in 
dirty, i.e. disordered, non-interacting electronic systems in two spatial dimensions, subject to static random impurities. 
Indeed, the aim of [5] was to study the integer quantum Hall plateau transition 1 , and to provide an alternative to the 
formulation given by Levine, Libby and Pruisken [7-1 1] in terms of a strongly coupled non-linear sigma model with 
a topological term. More recently it was recognized [12-15] that disordered superconductors (and other systems) 
provide an entire new arena capable of exhibiting (de-)localization transitions of a similar kind (albeit in entirely new 
universality classes). Since then, the study of (de-)localization transitions in non-interacting quantum systems in 2D 
has seen an immense surge of research activity. 

A general understanding of the strong-coupling (long-distance) behavior of 2D Kac-Moody current algebras 
perturbed by current-current interactions would be a valuable tool to describe a number of such transitions. Indeed, 
a good understanding of such perturbations has been achieved in a few cases [5,16-18]. However, in general, this is 
not the case, to date. 

An intriguing conjecture has recently been advanced by Gerganov, LeClair and Moriconi [19], who consider, as 
above, a general Kac-Moody current algebra conformal field theory with symmetry group Q at any level k, perturbed 
by right-left current bilinears. (The perturbations they discuss may also be anisotropic, or involve a supergroup, but 
this will not be important for the arguments presented in this article, which focuses entirely on the isotropic situation 
and bosonic groups.) Their paper builds on earlier work by Kutasov [20], who computed the renormalization group 
(RG) /^-function for the isotropic case 2 of a (bosonic) symmetry group Q to leading order in the large level k of the 
current algebra. The authors of [19] argued that Kutasov 's result be exact for any value of the level k, in a particular 
regularization scheme. Specifically, for the isotropic case, the authors of Ref. [19] conjecture that the exact (3-f unction 
for the coupling g be given by 3 

dg 1 C 2 g 2 

m := di= 2(TTW ' (L1) 

where C*2 is the eigenvalue of the quadratic Casimir operator in the adjoint representation of the symmetry group 
Q. Clearly, the notion of an exact [3-f unction is delicate due to its dependence on the regularization scheme. (The 
contributions to the /3-function beyond 2-loop order are scheme dependent.) Ref. [19] appears to be working in some 
scheme related to the left-right factorization of the underlying CFT. However, an explicit cut-off procedure, within 
which Eq. (1.1) is to be valid, is not specified in more detail in [19]. The authors indicate that certain checks to 3-loop 
order were performed. Checks beyond three loop order have never been performed, to our knowledge. The 3-loop 
/3-function within dimensional regularization has also been discussed in Refs. [21,22]. 

For the case where a symmetry Q = SU(2) is broken down to U(l) by a purely imaginary easy-axis anisotropy, 
and for level k = 1, the above conjecture (more precisely, an appropriate generalization) reproduces [23] known 
exact results [24-26]. The massless RG flow of the resulting non-unitary theory interpolates between ultraviolet 



'in the absence of long-range electron-electron interactions 
2 This theory is renormalizable with a single coupling constant. 

3 l := ln(a/L) is the RG-flow parameter, and a and L are the UV and IR cutoffs, respectively. 
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Figure 2.1: Group theoretical numbers for the classical groups, used in the main text. 



(UV) and massless infrared (IR) fixed points. (Both lie on the line of free scalar field theories with central charge 
c = 1 but with different compactification radii 4 ). In this case, the conjectured /3-function reproduces correctly the 
exactly known universal relationship between the exact scaling exponents (i.e. the slopes of the /3-function) at the IR 
and the UV fixed points. This is the only universal information contained in any /3-function describing this flow. 

On the other hand, the conjectured form of the /3-function, when appropriately generalized to supergroups, and to 
the anisotropic case, has recently been applied [28] to theories describing disordered systems, of the kind mentioned 
above. Here however, certain problems were encountered: Integration of the conjectured RG equations led to flows 
which reached a singularity after a finite scale transformation, which appears to be an unacceptable result. 

Motivated by these inconclusive results concerning the validity of the conjecture, we were led to check the conjec- 
tured form of the /3-function by explicit computation to high loop order. We consider the classical symmetry groups 
Q = SU(iV), SO(iV), and SP(iV), and the special case of isotropic current-current interactions at level k = 0. Our 
results are summarized in the following section. For all the classical groups, we find a contribution to the /3-function 
at 4-loop order which is incompatible with the conjecture in all possible regularization schemes. The discrepancy 
is caused by an extra logarithmic divergence in perturbation theory, proportional to an additional group theoretical 
invariant (besides the quadratic Casimir), which first appears at 4-loop order. This divergence is not accounted for by 
the conjectured form of the /3-function. Implications of our results, obtained for level k = 0, for the fc-dependence of 
the /3-function in any scheme are discussed in the conclusion, section 5. In this section, we also come back to, and 
comment on the special case of the anisotropic SU(2) model mentioned above. The reader who wishes to skip the 
technical details of our calculation, which are presented in sections 3 and 4, will find a self-contained exposition of 
our results in sections 2 and 5. 



In order to check the conjecture, we consider, as mentioned above, the specific case of an isotropic perturbation with 
symmetry group Q, and level k = 0. The conjectured /3-function (1.1) then becomes 



Here C 2 denotes the eigenvalue of the quadratic Casimir invariant in the adjoint representation of Q. For the classical 
groups, Q = SU(iV), SO(iV), and SP(iV), these are listed, in our normalizations, in figure 2.1. In sections 4.1 to 
4.7 we present an explicit perturbative calculation of the /3-function up to 4-loop order. This calculation proceeds in 
three steps: 

(i) use the current-algebra to calculate the diagrams, 

(ii) simplify the diagrams using elementary algebra, 

(iii) evaluate the integrals, which represent the ("Feynman") diagrams. 

After step (i), we encounter a great number of (rather complicated looking) diagrams. However after step (ii), we are 
left with only two classes of diagrams, 



4 Due to the non-unitarity of the theory, this is not in conflict with Zamolodchikov's c-Theorem [27]. 



2 Presentation of results 



M = ^C 2 g 2 . 



(2.1) 
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Figure 2.2: Chain-diagrams up to 4-loop order. The number of bubbles is the number of loops. 
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Figure 2.3: The group-theoretical invariant d,2 and its graphical representation. 



1. chain diagrams ("bubble" diagrams), 

2. non-chain diagrams. 

Chain-diagrams appear at loop-order 1, 2, 3, and 4. They have the form depicted in figure 2.2. Each bubble comes 
with a factor of g (the coupling-constant), with a group-theoretical factor of C 2 , and the whole chain with an (TV- 
independent) integral l n , at n-loop order (n bubbles in the chain). The integral Z n depends on the cutoffs L (infrared) 
and a (ultraviolet), and is a polynomial 5 of degree n in [ln(^)], plus terms which are finite for L/a — ► oo. (The 
fact that an n-loop integral is bounded by c [ln(^)] n with some constant c is a necessity to ensure renormalizability.) 
Only the leading term in Z n , with the highest power of ln(^), is universal. This applies e.g. to the diverging part 
of the 1-loop integral X\. Thus the contribution of the chain diagrams to the renormalization of the coupling g is, at 
n-loop order (up to a combinatorial factor) 



Non-chain diagrams first appear at 4-loop order. They are proportional to d 2 , which is an additional group-theoretical 
invariant (in the adjoint representation), independent of the quadratic Casimir C 2 . Its value for the classical groups 
is given on figure 2.1. This invariant can be constructed by drawing a cube, where one puts a factor of jf 3 on each 
corner, with one of its three indices on each adjoining edge, see figure 2.3. Finally, indices on the same edge are 
contracted. 

Each non-chain diagram has a global divergence proportional to [ln(^)] ("single log"), and subdivergences 
(higher powers of [ln(^)]). However, it turns out that to the order considered here, one can always group non- 
chain diagrams together into classes, such that each class has only a global divergence, but no subdivergence. This 
means that the integral (over positions) is proportional to ln(^) + finite, and that the prefactor in front of [ln(^)] is 
again universal. 

Let us recall here that a diagram proportional to [ln(^)], i.e. a "single log", gives a finite contribution to the 
/^-function. Using the above information, we thus find the following /3-function at 4-loop order 



Here, the numbers 03 and 04 depend on the regularization scheme (but not on N or the value of the cutoffs). In 
contrast, the remaining three terms are universal. The first term comes from the "single log" of the 1-loop chain 



5 See (A.7) for concreteness. 



g(C 2 g) n l n ■ 



(2.2) 



Pi9) = 9 \\C29 + a 2 {C 2 g) 2 + a 3 (C 2 gf + a 4 (C 2 g) 4 - d 2 ^(6 + tt v| + 0(g 6 ) . 



(2.3) 
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(giving the contribution \C2ti) and the last term from that of the non-chain diagrams (giving (6 + ir 2 )g 4 ). 

Furthermore, we have (12 = since we consider level k = 0; this is a consequence of the universality of the (3- 
function up to 2 loops, and will be checked for a specific scheme in appendix A.l. Note that to arrive at the result in 
(2.3), no specific form of the cutoff procedure has to be chosen. 

We are now in a position to answer the question, of whether the conjecture is compatible with our explicit 
calculation, in some given cut-off scheme. To see this, consider for example the group Q = SU(iV), where we have 
C2 = N, and d.2 = §iV 2 (see Fig. 2.1). Thus at 4-loop order, we have the following contributions 
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(2.4) 



Since C2 and d,2 contain all the dependence on N , and since 04 is independent of N, there is no possible choice of 
04, and thus no cutoff procedure, which cancels this term for all N. This proves, that the conjecture is incorrect for 
all possible cutoff procedures, for level k = 0. The same conclusion is arrived at for the other choices of groups. 

Let us finally give the result in a specific scheme, namely in the scheme in which an n-loop chain diagram is 
proportional to [ln(-^)] n , with no subleading term in ln(^): 

Pis) = \c 2 g 2 - d 2 ^(6 + vrV + 0(g e ) . (2.5) 

Indeed, this is the scheme used in the large-iV expansion of the Gross-Neveu model, where the /3-function becomes 
quadratic [29], to leading order in 1/N. Note that this is compatible with our result (2.5), since (upon rescaling g by 
N), the first (1-loop) term is order 0(1), whereas the second (4-loop) term is order 0(1/N 2 ). 

In conclusion, we have found that the conjecture is violated at 4-loop order, and at order 1/N 2 for SU(iV). For 
SO(iV) and SP(iV) corrections appear at order 1/N, as can be seen from the table in figure 2.1. 

Let us now outline the organization of the article; the reader wishing to skip the technical details of our paper 
can proceed directly to section 5: In section 3 we introduce the model, the current-algebra and basic notations. 
Our calculations are presented in section 4: We show in section 4.1 how the Kac-Moody current-algebra is used 
to successively eliminate interaction vertices from expectation values, and how this can be used to evaluate OPE- 
coefficients. This is a non-trivial task. Indeed, the raw result of this reduction procedure depends on the order of the 
successive reductions and is highly asymmetric, whereas the OPE-coefficient should be symmetric. To obtain a more 
symmetric result, the raw result can be simplified by using algebraic relations which we have baptized "magic rules", 
for their efficiency. This will explicitly be demonstrated in section 4.2 on the example of the 2-loop diagrams. In 
section 4.3 we proceed to 3-loop order, and show again how the initial highly asymmetrically looking OPE-coefficient 
is simplified. As in 2-loop order, all resulting diagrams are chain-diagrams which in a suitable scheme factorize, and 
thus do not give a new contribution to the /3-function. Proceeding to 4-loop order in section 4.4, one finds diagrams 
which the magic rules are no longer able to simplify to chain-diagrams. Due to the sheer number of initial diagrams, 
namely 576, this approach is not very illuminating. In section 4.5 we therefore pursue a different route: We first 
calculate OPE-coefficients for "adjoint" perturbations <I> a = f abc J b J c . We then show in section 4.6 how an n-loop 
OPE-coefficient can be expressed as a simple algebraic function times an OPE-coefficient at order n — 1, involving 
two adjoint perturbations. This allows us to identify at 4-loop order a combination of eight diagrams, which can not 
be factorized as chains (non-chain diagrams). Their contribution to the 4-loop /3-function is calculated analytically 
in appendix B. All these ingredients are collected in section 4.7, where we obtain the 4-loop /3-function. We have 
relegated some basic group-theoretical relations to appendices C.l, C.2 and C.3. Conclusions and further perspectives 
are offered in section 5. 



3 Model and method 

We study the Non-Abelian Thirring Model (NATM) in two dimensions. The model may be defined as a perturbation 
of a 2D conformal field theory, with action So, which is invariant under a symmetry group Q acting in the standard 
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way [1,30]. The chiral components of the conserved Noether currents, J a (z) and J a (z), depend (as indicated) only 
on z = x + iy and z = x — iy, and satisfy the defining operator product expansion (OPE) of the Affine Lie algebra 
(Kac-Moody algebra) at level k 

r( z ) = ^r + -ffJ c (o) + ... , J a (z)J b (o) = 1 ^r + ^fc b J c (o) + ... , (3.i) 

where ff 3 are the structure constants of Q. (Repeated indices are summed throughout this paper, unless stated 
otherwise.) The model we study is defined by the action 6 

S = S +goJ*(z,z), ^ : =y^, (3.2) 

where the perturbing operator 

&(z,z) = J a {z)J a {z) (3.3) 

is invariant under global transformations of the symmetry group Q. This theory is known to be renormalizable with a 
single coupling g. The conjectured form [19] of the /3-function for the renormalized coupling g is quoted in (1.1). 

We compute the /3-function explicitly to 4-loop order. To this end, consider the perturbative evaluation of the 
expectation value in the fully interacting theory of some quantity O, which may represent an operator, or a product 
of operators at different spatial positions, 

(01 »=^(°«- ! ' t,M ),' <«> 

Here, the expectation value (. . .) is taken in the unperturbed CFT with action Sq, normalized such that (1) = 1. 
Z(go) is the fully interacting partition function 6 . A cut-off (regularization) procedure, depending on short- and 
large-distance cut-offs a and L, is required to render all terms in this expansion finite, and is specified below. The 
renormalized coupling g (which depends on g$ and ^) is found by computing how the coefficient go of the first order 
term in the expansion of the exponential 

(••• l> + (...(-fl0jf +(•"(§/ J Z HZ1,Z1)HZ 2 ,Z2] 

+ (" ' / / / fl)*(z 2 , ^2)^(^3, Z 3 )j^+... (3.5) 

is modified by the higher order expansion terms. This modification is independent of the potential presence of 
any operator O in this expectation value, indicated by the ellipses 7 . The required calculation can be conveniently 
expressed in terms of multiple OPE coefficients of the perturbing operator $>(z,z), evaluated in the unperturbed 
theory. The product of (n + 1) such operators at different positions may be expanded into a "complete set" of 
operators & A sitting at the position of, say, the last operator. The expansion coefficients depend on the n relative 
coordinates, 

<f>(z 1 ,z 1 )...<$>(z n ,z n )<!>(z n+1 ,z n+1 ) = (3.6) 

= ^2 C ^[( Zl ~ Z n+l), (Zl ~ Zn+l); ■■■;(z n - Z n+1 ), (z n - Z n +l)] <& A {z n+ i , Z n +1 ) • 
A 

The non-vanishing expansion coefficients are exactly known in any CFT. In the present case they are especially 
simple, and can be obtained by successive use of the OPE of the currents, (3.1). In particular, the perturbing operator 



6 The partition function is Z (go) — f V [fields] exp(— S). 

The operator itself requires an analogous treatment, which, however, can be discussed independently; this will not be needed here. 
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<E> is the most relevant operator (besides the identity when k ^ 0) appearing amongst the § A ; all others are irrelevant. 
We find it convenient to denote the needed multiple OPE coefficient, where <E> A = by the symbol 



($(zi, Zl) . . . $0„, z n )$(z n+1 , z n+1 )\®(z n+ i, z n+1 )) 

= C$[(zi - z n+ i), (zi - z n +i); ...;(z n - z n+1 ), (z n - z n +i)] 



(3.7) 



The renormalization process is now easily understood by inserting (3.6) into (3.5). Explicitly, denote the relevant 
integrals over the multiple OPE coefficients by T n (for "Feynman"-diagram): 



Tn ■ = / {®{zi,Zl) ■ ■ ■ ®( z n+l, z n+l)\H z n+l, z n+l)) C(zi,Zi, Z n+1 ,Z n+1 ) 

J Zl,Z2,---Z n 



(3.8) 



These integrals are regularized by a cut-off prescription, which is achieved by inserting a cut-off function 

C[z\,z\, . . . , z n+ i,z n+ i) in the integral, as indicated. There are many possible choices. In this paper we choose a 

(circular) hard cut-off implemented by 



C(z 1 ,z 1 , ...,z n , Zn) := Y[ ( a < \ z i - Z j\ < L ) 

i+3 



(3.9) 



where is the usual step function. This cut-off procedure restricts the distances between any pair of integration 
variables to lie between the short- and the long-distance cut-offs a and L. All integrals T n are thus finite functions of 
■g. As usual, inserting (3.6) in (3.5), and using (3.8) gives: 



( ...i>„ + (...(-»/.(x,fl)) +(...(!«/.(«,. 



+ 



-So 3 



3! 



+ 



= (...l)o+ ... 



' / Hz,. 

J z 



+ ... = 



-gf z $(z,z) 



Following standard reasoning we have re-exponentiated in the last line. One can now read off the renormalized 
coupling: 



(3.10) 



The /3-function is obtained as the change of g in response to changing a (or 1/L), while keeping the bare coupling go 
fixed: 



9 



90 



( 90 > z) 



(3.11) 



In the remaining sections of the paper we will obtain the integrals T\,...,T±. This gives us the result for the 4-loop 
/^-function written in (2.3) above. 



4 Calculation 



In this section we present in detail the evaluation of the integrals T n defined in (3.8) ("Feynman"-diagrams), needed to 
obtain the /3-function, as explained in Section 3. The core of this calculation consists in obtaining the OPE coefficients 
defined in (3.7), by repeated use of the current algebra OPE (3.1). We start with the simplest case, i.e. with the 1-loop 
integral T\, and proceed successively to the more involved cases, up to 4-loop order. 
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Figure 4.1: Labeling of the points in figures 4.2, 4.3, 4.4, and diagrams in the main text. 



4.1 1 -loop order 

At 1-loop order, we need the OPE-coefficient 

($(z,z)&(w,w)\$(w,w)) . 
To evaluate it, we have to eliminate &(z, z) from $(z, f w), using (3.1). This is done as follows 8 

$(z,z)$(w,w) = J a {z)J a {z)J b (w)J b {w) 

— ffffj c (w)J d (w)-^—, = C 2 J c (w)J c (wh 



\w — z\ 



\w — z 



2 ' 



(4.1) 



(4.2) 



We have used that f abc f abd = C2<5 cd , with the second Casimir Ci- (This and more group theoretical relations are 
derived in appendix C.l.) We denote (4.1) in short by 



z)$(w,w)\&(w,w)) = 



Co 



w — z 



12 • 



(4.3) 



The arrows show the direction in which the elimination has been made. This defines the sign. To be specific, an arrow 
from z to w represents l/(z — w). A dashed such arrow represents l/(z — w). When a solid and a dashed arrow 
(with the same direction) connect the same two points, one can drop the arrows for simplicity of notation; seeing a 
solid and a dashed line thus means that when adding the arrows, both arrows are pointing in the same direction. 
The OPE-coefficient (4.3) yields the 1-loop diagram T\ 



F 1 = l<^=l ] ^ w e(a<\ W -z\<L) = C 2 ln^). 



(4.4) 



4.2 2-loop order and the magic rule 

At 2-loop order, we have three $'s. Denoting $j := J a (zi)J a (zi), we need to calculate ($i$ 2 ^3 |$3)- Straight- 
forward use of the current-algebra (3.1) with k = 0, eliminating the currents one by one, starting with point 1, and 
continuing with point 2, yields 



($i$ 2 $ 3 |$ 3 ) = ifl 



|Z12| 2 |Z23| 2 



+ 



■^13 1 2 1 ^23 1 2 



R23 ^13^12 



1^231^12^13 



(4.5) 



where we have abbreviated z. 



Zi — Zj. Here, and throughout this article, we use the labeling of points as 
indicated in figure 4.1. The result is graphically presented in figure 4.2 (top). Equation (4.5) apparently contains a 
new diagram, which renders the OPE-coefficient asymmetric upon exchange of point 1 with point 2, or of point 1 
with point 3. However, there is a simple algebraic identity, the "magic" rule for the real part 3? of 



zw ' 



ft 



" 1 " 


= ft 


zw 




zw 
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zw 


\z 2 \w 2 




_ \z\ 2 \w 2 


~ 2 



1 

w 



+ 



w 



\z — w\ 2 
\z\ 2 \w\ 2 



(4.6) 



Recall that summation over repeated indices is implied, and that we work at level k = 0. 
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Figure 4.2: 2-loop diagrams after reducing the structure-constants to numbers. To be multiplied by \Cl- The numbers given are the weight. 
The first line is the raw result, as obtained by using the reduction algorithm. Arrows indicate the direction of the reduction. The second line 
after using magic relations. 



The most useful application is in the presence of an additional factor l/\w — z\ 2 , which cancels the numerator in the 
last term. This leads to the decomposition of the new diagrams in (4.6) into chain-diagrams [drawn below rotated by 
—120° as compared to figure 4.2 (top)] 



(4.7) 






The OPE-coefficient (4.5) simplifies to 

($i$ 2 $3|$3) = \C\ 



Ul2| 2 Ul3| 2 



+ 



+ 



^13| 2 k23| 2 



(4.8) 



which is manifestly symmetric, as it should be. The resulting expression for the OPE coefficient in (4.8) is graphically 
represented in figure 4.2 (bottom). One sees that after using the "magic" rule, the OPE coefficient, and hence the 
integral JF 2 , can be written in terms of chain diagrams. This suggests, that the corresponding diagrams (i.e. the 
"Feynman" integral Ti) factorize, are of order ln(-^) 2 without a pure ln(-) and thus give no contribution to the 
/3-function at 2-loop order. This is indeed correct, as checked in appendix A. 1 for the cutoff -procedure introduced in 
section 3. 

For the model at hand, the cut-off procedure is subtle. The reason is that one cannot put a cut-off on the lines, 
as would be most convenient to immediately prove factorization of chain-diagrams: In constructing the diagram, we 
have used magic rules to move around the lines, and if we leave behind a cut-off function, then the resulting diagram 
will not be totally symmetric, as it should and as it is in our construction. The only way out of the above dilemma, 
is to put cut-offs between any pair of points, regardless of whether the two points are connected with a line or not 
[compare (3.8)]. However then the factorization is no longer a trivial statement, and has to be checked. This has been 
done for the 2-loop chains in appendix A.l. As we have argued in section 2, this is not essential for our arguments, 
and the conclusions remain valid in any scheme. Let us however mention, that in order to recover the large- N limit 
of SU(iV), factorization is needed, and is sufficient to uniquely fix the RG-procedure up to 4-loop order; but not 
necessarily beyond. 



4.3 3-loop order 

At 3-loop order, 36 diagrams appear, presented on top of figure 4.3. These diagrams all contain six structure-constants, 
and have two free indices a and b, which are contracted with the remaining J a J b . Since the only invariant object 
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magic rule 




Figure 4.3: 3-loop diagrams after reducing the structure-constants to numbers. To be multiplied by \C\. The numbers given are the weight. 
The first 4 lines are the raw result, as obtained by using the reduction algorithm. Note that 4 diagrams have weight 0. Arrows indicate the 
direction of the reduction. The last two lines after using magic relations, dropping the redundant arrows. 

with two indices is 5 ab , one can contract the last lines to obtain the algebraic factor; the final result has of course to 
be divided by the dimension of the adjoint representation. One can then convince oneself by drawing pictures, that 
all objects which can be constructed, contain at least one loop made out of two or three vertices 9 , i.e. objects of the 
form 




= -C 2 S 



ad 



jaed jbge jcdg 



G 



2 rabc 



(4.9) 



(4.10) 



which we give together with a group-theoretical identity (derived in appendix C.l), which is sufficient to reduce the 
number of / in any given diagram. Repeatedly using (4.9) and (4.10) thus allows to eliminate all /. This procedure 



The simplest object without such a loop would be a cube, which indeed appears at 4-loop order, see figure 2.3. 
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Figure 4.4: 4-loop diagrams not proportional to C| . 



is performed using a computer, and the reader would have a hard time verifying it by hand. We have thus shown that 
all 3-loop diagrams are proportional to Cf, thus no additional group theory invariants, besides the second Casimir, 
appear at this order. 

The diagrams are given with their combinatorial factor on top of figure 4.3, to be multiplied by ^Cf. Applying 
magic rules leads to chain-diagrams, presented graphically at the bottom of figure 4.3. Algebraically, the result is 



C 3 

($l$2<£3$4|$l) = -1 



1 



+ 



1 



+ 



1 



|^12| 2 |^14| 2 |^23| 2 |Z13| 2 |<Z14| 2 |Z23| 2 1^12 1 2 1 Zl3 ? 1^24 1 2 
1 1 1 



+ 



+ 



+ 



kl3| 2 kl4| 2 k24| 2 kl3| 2 k23| 2 N24| 2 1^14 1 2 1^23 1 2 1^24 1 2 

1 1 1 

1 ~\ Tol To 1 Tol m To I Tol fol To" 

F12TF23I F34T 



Ul2| 2 Ul3| 2 k34| 2 



kl2| 2 kl4| 2 |^34| 2 



+ 



+ 



+ 



^14| 2 k 23 | 2 k34| 2 M 2 k24| 2 N 4 | 2 kl 3 | 2 k24| 2 k34| 2 



(4.11) 
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4.4 4-loop order, direct approach 



Let us now continue to 4-loop order. After using the current- algebra, there are 576 diagrams, which again we generate 
computer-algebraically. The group theoretical factors appearing with these diagrams are much more involved. An 
example is a cube, where each corner represents a structure-factor f abc and each link identifies a pair of common 
indices between two /'s. This is drawn on figure 2.3 and detailed in appendix C.l. It is at this loop-order that 
an additional group theoretical invariant besides the quadratic Casimir arises. After reducing the algebra, one finds 
that 380 terms are proportional to C|. Using magic rules, these diagrams can be reduced to 60 chains; these are 
in fact all the chains which can be drawn through 5 points 10 . Each chain comes with a weight of |C|. For the 
remaining diagrams not proportional to C| , presented in figure 4.4, our reduction-algorithm based on magic rules is 
incapable of further simplifying it. In section 4.6 we will present a simple calculation, reducing the task to calculating 
a combination of eight diagrams. To this aim, we need correlation functions involving operators which we call 
"adjoint" perturbations, defined below. 



4.5 OPE for adjoint perturbations 

Define the "adjoint" perturbation at position (zi, Zj) as 

^{ Zi ,Zi) = ^1 ■= f abc J b (z t )J c (z t ) . 



(4.12) 



We now apply the same procedure as in the previous sections: Eliminate currents one by one using i) the current- 
algebra, ii) evaluation of the group theoretical factors, and iii) simplifications with the magic rule. After some lengthy 
calculations (done again computer-algebraically), we find up to 3-loop order: 



($«$»|$ 2 ) 
($«$ 2 $«|$ 3 ) 

($«$«$3$ 4 |$ 4 ) 



c 2 



1 



2 M 



kl2| 2 kl3| 2 



+ 



Zl^\ 2 \z2i\ 2 



kl2| 2 kl4| 2 k23| 2 



+ 



^12| 2 |^13| 2 k24| 2 



+ 



^12 | 2 |^13 I 2 |^34| 2 



+ 



^12| 2 Ul4| 2 |«34| 2 



+ 



^12| 2 k23| 2 |«34| 2 



+ 



Z12\ 2 \z2a\ 2 \zm\ 2 



+d 2 





+ 




(4.13) 
(4.14) 



(4.15) 



In 2 



com- 



The additional group-theoretical invariant cfe is defined 11 in appendix C.l. For SU(iV) this reads d 2 
pared to the leading term C| = N 2 . The results for SO(iV) and SP(iV) are listed in figure 2.1, see also appendix C.3. 
For these groups, C| ~ iV 4 and again c?2 is subdominant, with cfo ~ N 3 . 



10 In general, there is a total of |n ■ (n — 1)! chains that can be drawn through n points. 

"Constructing a symmetrized tensor d abcd out of the trace of 4 structure-constants /" b , d,2 is the non-dominant contribution (in TV) of the 
square of d abcd [as defined in (C.13)]. 
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Figure 4.5: The combination of 8 diagrams contributing at 4-loop order. Note the different labeling of points given in the inset, as compared 
to the labels of figure 4.1, used in figure 4.4. 

4.6 4-loop order simplified 

We have seen in section 4.4 that a direct 4-loop calculation is quite cumbersome. Instead, we use here a different ap- 
proach, inspired by the original work by Kutasov [20] . We start by eliminating $ n from the multiple OPE coefficient 
(3.7). Let us first give the result and then explain how we have obtained it: 

($ n $l$2---$n-l|$n-l) = 




(4.16) 



We have eliminated all currents at point n. Using the current-algebra (3.1) again with k = 0, there is a contribution 
from each pair of points {i, j} with i,j / n. The first line of (4.16) contains the contributions with i / j, for which 
we have listed below the corresponding current-algebra identities in (4.17) and (4.18). The last line of (4.16) is the 
case i = j, and is obtained by using the current-algebra both for the holomorphic and antiholomorphic current, as 
given in (4.19) below. 

— (4.17) 

Z-n Zi 

1 _ $° (4.18) 

Z n Zj 

= . C * . $i ■ (4.19) 

\Zn Z%\ 

Note that eliminating point n (instead of point 1 as we were used to do) is for later calculational (and representational) 
convenience only. 

We now turn to the 4-loop calculation, i.e. set n = 5. One can check that starting from (4.16), using (4.11) 
and (4.15), one reconstructs all the 60 chains connecting 5 points, as found in section 4.4. The remaining terms are 
obtained from the first term in (4.16) times the term proportional to c?2 in (4.15). There are ( n 1 K" 2 ) = g suc h terms, 
each being a combination of 8 diagrams depicted in figure 4.5. Since each of the 6 terms gives the same contribution 



J a (z n ) J*(^)J*(Zi) 



cabc 



r(z n ) J b {z 3 )j\ Zj ) 



Zyi z% 
cabc 



J a (z n )J a (z n ) J b (z t )j\ Zl ) 



f° 



z. 

abc xabd 



■JTz 3 )J\z,) = 



r c f 



J C {Zi)J d {Zi 
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upon integration, we only have to calculate the integral over one of them. Analytically, this is most easily written as 
(we have chosen i, j to be the pair 1, 2 and the starting point is 5) 

l! = /(—-—) (— )t J t4— + — ) 

J \ZUZ13 ZuZ 13 J \Z24Z23 ^24^23/ \Zu\ V z 15^25 Z 15 Z 2 5 J 

xC(zi,zi, . . . ,z 5 ,z 5 ) , (4.20) 

where the integral is over all but one point, and the cutoff-function C(. . .) was introduced in (3.9). This integral is 
evaluated in appendix B to be 

1= - — (Q + 7T 2)l n (L\ . (4.21) 



12 v \a 



4.7 The /^-function up to 4-loop order 



Now we are ready to put everything together to obtain the /3-function. In a scheme in which the chains factorize, we 
obtain by collecting the results (4.4), (4.8), (4.11), and the paragraph below (4.16), and upon use of (3.10): 

9 = So jl - 7y [soC 2 ln(£)] + I| [ 90 C 2 ln(^] 2 - IH [ 5o C 2 ln(§)] 3 

+ [soC 2 ln(§ ) ] 4 + 6^(6 + . 2 ) ( i 2 ^ln(|)) |. (4.22) 

Note that the factors 1/n! are from the exponential, then for the chains the next factor is the number of chains times 
their dependence on go, times the group-theoretical factor 12 C% /2 n ~ 1 . The last term (which comes from the non- 
chain diagrams) has a factor of 6 from combinatorics as discussed in the previous section and the minus signs from 
the integral (4.21) and from (4.16) cancel. 

Inserting (4.22) into (3.11) leads to the 4-loop /3-function in terms of the renormalized coupling g: 

0(9) = \C29 2 - d 2 ^(6 + vrV + 0{g & ) , (4.23) 

with d 2 = fiV 2 for SU(JV), d 2 = 24 - fN + 6N 2 - fiV 3 for SO(N) and d 2 = | + § N + fiV 2 + -^N 3 for 
SP(iV), as calculated in appendices C.2 and C.3. 

In schemes in which the chains do not factorize, there are additional terms, see (2.3) and the discussion below 
that equation. 

Some comments on the procedure are in order. Readers used to the Wilson-scheme, will recover that procedure 
by studying the change of g in (4.22) under an infinitesimal change of a, corresponding to the integration over an 
infinitesimal shell from a to a + 5a. The only difference is that this is a shell in position space, and not in momentum 
space. 

Second, to our knowledge this is the first 4-loop calculation with a hard cutoff, or equivalently the first 4-loop 
calculation in a Wilson scheme. 

5 Conclusion and further perspectives 

In this article we have performed an explicit perturbative calculation of the /3-function for the non-abelian Thirring 
model at k = up to 4-loop order. We have found that the conjectured form of the /3-function [19], Eq.(l.l), is 



See (4.4), (4.8), (4.11), and the paragraph below (4.16). 
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incompatible with our result in all regularization schemes. The discrepancy arises from an extra logarithmic diver- 
gence, which appears first at 4-loop order, and which is proportional to a higher group-theoretical invariant (evaluated 
in the adjoint representation of the symmetry group) which is different from the quadratic Casimir invariant. This 
divergence is not accounted for by the conjectured /3-function. 

It is worth pointing out that our explicit 4-loop result at level k = does not only rule out the particular conjec- 
tured form of the (isotropic) /3-function Eq. (1.1), but a more general class of conjectures for the /3-function. This 
way of presenting our 4-loop result emphasizes the dependence on the level k, whereas in section 2 only the special 
case k = was discussed. Such forms which we can rule out arise by attempting to "scale" with the level k. Specif- 
ically, for any one of the classical groups Q = SU(iV), SO(iV) and SP(iV), the /3-function (of the isotropic theory) 
will in general be a function of three variables, the coupling constant g, the level k, as well as N, or equivalently 
C 2 = C2(N), the second Casimir invariant in the adjoint representation: 

ijL = P(g,k,C 2 ). (5.1) 

It was argued in [23] that by rescaling the Kac-Moody currents J a — > \fkJ a (as suggested by the large-/c calculations 
done in [20]), the (isotropic) /3-function should satisfy (perhaps in a suitable scheme) the "scaling form" 

l3(g,k,C 2 )i^F(^-,kg] =g H (gC 2 ,^-\ (5.2) 

The conjectured /3-function of Ref. [19], i.e. Eq.(l.l), is a special case of this. The second equation above gives an 
equivalent way of writing the scaling form, useful when considering the limit k — > for g = fixed, whereas the form 
in the first equation is useful in the large k limit where 1/k —> for kg = fixed. Since we know that the perturbative 
/3-function must have a finite limit as k —> 0, the second equation in (5.2), when specialized to k = 0, leads to a 
form of the /^-function, whose ^-dependence is only through the combination gC 2 (apart from an overall factor of g). 
Comparison with Eq. (2.3) shows that this is incompatible with the explicit 4-loop result that we have found in any 
possible scheme. Hence, our result implies that the /3-function must have an explicit dependence on the level k, and 
that the latter can in no scheme be "scaled out" in the way indicated in (5.2). 

Finally, a point which deserves clarification is why in the case where a symmetry Q = SU(2) is broken down to 
U(l) by a purely imaginary easy-axis anisotropy, and for level k = 1, the conjecture reproduces [23] known exact 
results [24-26]. (This was mentioned in the introduction, section 1.) Indeed, Ref. [23] proposes that this agreement 
should provide a strong check of the conjecture. Here we would like to point out, however, that this agreement is 
not surprising, because this theory is very special, and highly constraint by its underlying hidden quantum group 
symmetry (or, "fractional supersymmetry") [31], present for all values of the level k. The symmetry imposes strong 
constraints on the ^-dependence of the relationship between the slopes of the /^-function, i.e. the RG eigenvalues y 
of the perturbation, at the UV and IR fixed points. As a consequence of the symmetry, this relationship is [32] : 

' + r^— = 1 . (5.3) 



k ym k yuv 

Note that the same result can also be obtained by using the conjectured (3-t unction of [19], see [23]. Now since, 
following Kutasov [20], the conjecture is the leading term in a 1/ fc-expansion, it should yield a relation between ym 
and y\jY, which is valid at leading order in 1/k, but will not contain information about corrections to this of order 
1/ A; 2 or higher. However, due to the fractional supersymmetry, the exact relation between yiR and yuv has no such 
higher order corrections at all, as can be seen from Eq. (5.3). Thus, due to the constraints imposed by this symmetry, 
the leading order term in 1/k happens to give already the whole, i.e. the exact result. This explains why the conjecture 
reproduces the exact result even for level k = 1. 

We end our discussion by noting that it would be interesting to obtain, generalizing Kutasov 's work, who com- 
puted (as mentioned) the /^-function of the Non-Abelian Thirring model to first order in 1/k in the large-/c expansion, 
higher order terms in this large-/c expansion. These will not in general be absent, as our work presented in this paper 
shows. Work along these lines is in progress. 
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Appendices 



A Factorization of chain-diagrams 

In this appendix, we show how chain-diagrams factor, restricting ourselves to 2-loop order. This is done in appendix 
A.l. As a tool, we need the "conformal mapping technique", which was introduced in [33,34], reviewed in [35], and 
which we present here for completeness, and since in contrast to the cited references we here work exactly at the 
critical dimension, where we need both an ultraviolet and an infrared cutoff. 



A.1 Factorization of chain-diagrams at 2-loop order 



At 2-loop order, everything can with the help of the magic relation be reduced to the bubble-chain. The subtracted 
2-loop diagram, i.e. the 2-loop diagram minus the square of the 1-loop diagram is (we denote by S this subtraction- 
operator, which also contains the integration and the cut-off functions) 




1 



[0(a < \z\, \w\, \z — w\ < L) — 0(a < \z\, \w\ < L)\ 



(A.l) 



The first term on the r.h.s. represents the 2-loop integral, the second term the subtracted 1-loop integrals (where 
integration over w and z factorizes). Applying — a to the above gives 

9 s 

oa 




a L,w M 2 IH 



0(a = \z\ < \w\, \z — w\ < L) 

+ 0(a = \w\ < \z\, \z — w\ < L) + 0(a = \z — w\ < \w\, \z\ < L) 
- 6(a = \z\ < \w\ < L) — 9(a = \w\ < \z\ < L) . 



(A.l) 



Using the conformal mapping technique of [33-35], which is summarized in appendix A.2, all terms can be mapped 
onto \z\ = a; with the result (we have used that a 2 /\z\ 2 = 1): 




[9 ON = a, \w\, \w - z\) < §) - (Sg(a, \w\) < §)] 



The function i (ai, . . . , a n ) is defined as 



max 



_ max(ai, . . . ,a n ) 

nun v a i> • • • i a n) '■— •„/ n \ 

min(ai, . . . , a n ) 

The above is a function of L/a, and can be bounded for L/a > c\ by 



9 c 

oa 




a 

<C2 I 



(A.3) 



(A.4) 



(A.5) 



with C2 = 02(01). Taking c\ > 3 allows the bound C2 = 2. The important thing is that the integral is not diverging: 
this means we have subtracted the right 1-loop counter-term. Moreover, the limit of large L/a can be taken; since it 
is zero, there is no single ln-contribution in the 2-loop integral. We can denote symbolically the result as 

2 

(A.6) 
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A.2 "Conformal mapping" 

As a tool to prove factorization of chains (see appendix A.l), we need the "conformal mapping" technique, which 
was introduced in [33,34], reviewed in [35], and which we present here for completeness, and since in contrast to the 
cited references we here work exactly at the critical dimension, where we need both an ultraviolet and an infrared 
cutoff. 

Note that a general iV-loop integral In will behave as 

l N (a, L) = a + 01 In § + a 2 (in |) 2 + . . . + a N (in %) N , (A.7) 
where we dropped terms which vanish in the limit of L/a — > oo. Deriving w.r.t. a leads to 

-a—l n (a, L) = ai + 2a 2 In § + . . . + Na N (in . (A.8) 

On the level of the integral, this operation amounts to fixing the smallest distance to be a. Due to our normalizations, 
this is equivalent to fixing the both endpoints of this smallest distance. The integration over the remaining points has 
then to be done. 

We now state a very important theorem for the integral over a function / at order N — 1 loops : If / (z± , z\ , . . . , , z~n ) 
is a homogeneous function of dimension —2(N — 1) (z and z have dimension 1), then the integral over z\, . . . , zjv-i 
(the relative coordinates between points) 

X N (a,L):= f(z 1 ,zi,...,z N ,z N )C(zi,z l ,...,z N ,z N ) (A.9) 

J Z\,...,ZN-1 

has dimension 0. Consider a sector S (ordering of the distances). Bei Q := \z\— Zj\, with 1 < a < m := N{N—\)/2. 
Then S := {z±, . . . , zn}, s.t. x\ < x 2 < . . . < x m . (Actually, we have chosen the labeling of the distances x a to 
account for the ordering. This is not always the most practical thing to do.) Also define the characteristic function 
Xs(xi, . . . , x m ) of a sector S as being 1 if all distances satisfy the inequalities of the sector and otherwise. The 
a-derivative of the integral restricted to the sector S is 

J s :=-a—l N {a,L) = J f( Zl , . . . , z N )\ xi= @{x m < L) X s(xi, . . . , x m ) . (A. 10) 

The conformal mapping theorem [33-35], whose proof we reproduce below for completeness, now states that if the 
integral (A. 10) is Riemann-integrable everywhere, then 

J s = J f(z!,. . . ,z N )\ x _ a Q(x m /x 1 < L/a)xs(xi,...,x m ) . (A.ll) 

In words: The above integral can be evaluated by fixing any of the distances to be a (or 1 equivalently). The constraint 
on the smallest and largest distances is captured by the condition that the ratio of largest to smallest distance is 
bounded by L/a, as it is in the original integral, which is thus just a special case of the expression (A.ll). 

Proof: 



(A.12) 



First of all, since x\ = a, and introducing a (5-function to enforce it, J s becomes 

J S = j f(zi, ■ ■ .,z N )5(x 1 - a)Q(x m /x 1 < L/a)xs(xi, ■ ■ -,x m ) . 
We now aim at integrating over distances x±, . . . , x m instead of coordinates with an arbitrary function g 

/dV...dV_ 19 (, 1 ,...,x m ) = /dx 1 ...d, ra ,(x„...,x ra Mx„...,x m ). (A.,3) 
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The measure is easily constructed as 

fi(xi, . . .,x rn ) = j d 2 zi . . .d 2 z N _i 5{xi - \zi - z 2 \) . . . 5(x m - \z N -i - z N \) , (A. 14) 

where the (5-distributions enforce the Xj's to be the distances between the Zj's. 

We now want to map onto x\ = a. To achieve this, we can always do the integration over x\ last. This gives for 

J S 

J s = J dxi J dxi . . . dx/_!dx/ + i . . . dx m fx(xi, x m ) 5(x 1 - a) 

x f(xi, . . .,x m ) @(x m /x 1 < L/a) xs(xi, ■ ■ -,x m ) . (A.15) 
We now make a change of variables. For all % but I, set 

Xi := XiXi/a . (A. 16) 

We also define xi := a, and introduce this into (A.15) as 1 = J* dxi 5{xi — a): 



J S = J dxi J dxi . . . dx m fj,(xi, x m ) 5{xi - a) 



x f(xi, x m ) Q(x m /xi < L/a) xs(xi, x m ) 

x 8{x\xi — a) — . (A.17) 

Note that the factor of ^ consists of (H^^ -1 )/ 2-1 f rom the terms dxi but dxf, a factor of (Si)^ -1 ^ 2- "?) f rom 
the measure; and a factor of (^) 2iyN ^ from /. Using that 

/dxi 5(x\Xi — a)— = 1 , (A.18) 
xi 

we obtain 

J S = j dxi . . . dx rn fi(x 1 , x rn ) 5(xi - a) 

x f(xi, . . . ,x m )Q{x m /xi < L/a)xs(xi, ■ ■ .,x m ) ■ (A. 19) 

Dropping the tildes, this is nothing but (A.15) with x\ replaced by x\ which completes the proof. 

B The 4-loop integral 

In this Appendix we evaluate analytically the integral (4.20) needed in Section (4.6) to obtain the universal part of 
the 4-loop contribution to the /^-function [i.e. the last term in (2.3)], with the result quoted in (4.21). The integrand of 
the integral (4.20) in question is 

M~(^-^-)(-^-^-)'(-± + J-). (B.1, 

\ZUZ 13 ZuZtfJ \Z24Z23 ^24^23/ \Zu\ V z 15^25 ^15^25/ 

Graphically, this is depicted in figure 4.5. We observe that we can make the following simplification (due to the 
"second magic rule"): 

1 1 urn — wu „. wxu n .\w — u\h ^ N 
— - — = - - = 2z , m , 2 = 2i \ , (B.2) 

wu wu wwuu \wr\ur \w\ z \u\ z 



18 



where h is the height of the triangle spanned by w and u; if the angle is larger than it, then h is negative. Graphically, 
this can be visualized as 



(B.3) 




I w - u I 



1 

NIP 



Note that the first two factors of the integrand M. both contribute a term | Z34 1 , thus canceling the third term 
This allows us to see that the integral has no subdivergences; it will contain only a "global divergence", i.e. it will 
be proportional to a single power of ln(L/a) (L and a are the IR and UV cutoffs, respectively). We now proceed to 
check this by explicit calculation and to compute the precise coefficient of the single logarithmic divergence. Let us 
now introduce distances as depicted in the figure: 




(B.4) 



Here all distances are measured from except for x' and z which are measured from their intersection point. In these 
conventions, x' and — z in the figure are negative. The integrand can then be written as 13 



M = {2if 



—z 



(y 2 + x 2 )(y 2 + (x> + 6) 2 ) (z 2 + x' 2 )(z 2 + {x + b) 2 ) V*i5*25 z 15 z 25 



+ 



(B.5) 



where all variables are to be integrated over. Let us first do the integrals over x, x', b, z and z^ = ^t(z^) + ^(25), i.e. 
all distances except for y, which is kept fixed and positive. At the end, we integrate over the vector y, both over its 
magnitude and direction. (This fixes the coordinate system.) We note that choosing a <^ \y\ <^ L, boundary terms 
can be neglected, since the integrals do not contain subdivergences, neither in the UV nor in the IR. 
Doing first the integral over point 5, we obtain using (D.l) from appendix D 



/ 



d 2 z 5 



1 



+ 



1 



^15^25 ^15^25 



^ = -2vrln(|zi 2 | 2 ) + const. = -2vr In [(y + zf + b 2 ] + const. , (B.6) 



where the constant depends on the IR-cutoff L. However one easily sees that it drops from the above calculation due 
to the asymmetry of z —> —z of the remaining terms in (B.5). We finally have to integrate: 



-8vr 



y 



{y 2 + x 2 ){y 2 + [x> + 6) 2 ) ( z 2 + x' 2 ){z 2 + (x + b) 2 ) 



In [{y + z) 2 + b 2 ] . 



(B.7) 



There are four complex integration variables, equivalent to eight real integration variables. We make use of this equivalence whenever 
convenient. 
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The simplest integrals are those over x and x' . We use 

1 1 



f 

j — i 



dx ■ 



ir(\y\ + \z\) 



x 2 + y 2 (b + x) 2 + z 2 \yz\(b 2 + (|y| + \z\) 2 ) ' 
which can be done by residue-calculus. Integrating (B.7) over x and x' thus gives: 



-8tH 



\y\ + \A 



■ln[(y + z) 2 + 6 2 ] . 



(B.8) 



(B.9) 



yz(b 2 + (\y\ + \z\) 2 ) 2 

To continue, we recall that by construction y (which is the module of a vector) is positive. (B.9) can thus be written 
as the integral over positive z only 

(y + z) 2 



-8tt j 



(In [{y + zf + b 2 ] - In [ { y- z ) 2 + h 2 ]) 



yz(b 2 + (y + z) 2 ) 2 

The easiest integral to do is that over b, which nevertheless is a little bit tricky. We need 



m(6 2 + d 2 ) 
(b 2 + s- 



2\2 



db 



— IT 



s 2 (\d\ + \s\ 



+ 



7rln(|d| + |s| 



(B.10) 



(B.ll) 



which can be verified with the help of the residue-theorem. To do so, one splits the ln(6 2 + d 2 ) = ln(6 + i\d\) + 
ln(6 — i\d\) which both have branch-cuts. But the integral can be closed either in the upper or lower domain, and we 
close it in the domain where there is no branch-cut. This leaves us with 

(y + z? 
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Scaling out y, and splitting the integral into domains where the absolute values have a definite sign gives: 
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The final integral over y contains the integral over the modulus of y and its direction, which contributes a factor of 
2tt: 



I d y 2 ^ ["^ 4 (6 + ^ 2 )^] =-^ 5 (6 + vr 2 ) ln(|) . 



(B.14) 



To conform to the normalizations used in the main text, see Eq. (3.2), this still has to be divided by (2ir) A , yielding the 
final result [with the integral running over all but one of the points, and normalizations according to equation (3.2)] 



J:= J M = -^(6 + vr 2 ) In finite 



(B.15) 



We have indicated an additional finite term in the result, which depends on the specific regularization prescription, 
and which is either a constant or decays to in the limit of L/a — > oo. 
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C Some remarks on group theory 



In this appendix, we collect a number of useful group-theoretical identities, first in appendix C.l for a general Lie- 
group Q, then in appendix C.2 for SU(JV), and finally in appendix C.3 for the other classical groups, SO(iV) and 
SP(iV). 

C.1 Group theoretical invariants 

In this appendix we discuss the additional group theoretical invariant, referred to in the main text. Since we are using 
the current algebra, only the adjoint representation of the symmetry group Q appears in our calculations. Therefore, all 
group-theoretical invariants that can possibly appear, can all be constructed out the structure constants. The simplest 
such invariant is of course the eigenvalue of the quadratic Casimir invariant C2 in the adjoint representation, which 
is of second order in structure constants f ahc . Here we consider invariants which are of higher order in the structure 
constants. 

Notation: The zero modes j a := = j(dz/2ni) J a (z) of the Kac-Moody currents [1] are the generators of 
the Lie-group Q, satisfying the commutation relations 

= fc ab f , (C.l) 
which are represented in the adjoint representation by matrices 

(T a )c b := f c ab (C.2) 

We work with antihermitean generators j a , so that the structure constants f c ab are real. 
The (/-invariant Killing form i] ab , and its inverse r} oc , defined by 

rf b := ^ tr {r a T b ) , V ab Vbc = 5 a c (C.3) 

may be used to raise and lower adjoint indices a,b,.... Here M is a suitable normalization constant. Then, (C.l) 
and (C.3) imply that the structure constants f cab = f abc are totally antisymmetric. Throughout this subsection, we 
choose a basis of the Lie algebra for which rf b = 5 ab . Hence, no distinction between upper and lower adjoint indices 
has to be made. [The matrices (T a ) b in (C.2) are then antihermitean.] 

We now proceed to discuss various group-theoretical invariants, needed in the main text, which can be constructed 
out of products of structure constants. Our discussion is organized according to the number of factors f abc appearing. 

Quadratic Casimir. The eigenvalue C2 of the quadratic Casimir invariant in the adjoint representation 14 , 

( T a T a^d = jab ^ad = gd ^ fabc jabd = ^ g cd (Q4) 

is of 2nd order in the structure constants. Eq. (C.4) is graphically depicted in (4.9). 

"Triangle Rule": The Jacobi-identity implies the following relation for the structure constants: 

fe ad fd bC + fe bd fd Ca + fe Cd fd ab = . (C.5) 

which is just (C.l): 

([T\T b \Y e = f d ab (t^" (C.6) 



As usual, all repeated indices are summed. 



J ,3 
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Multiplying (C.5) with f\ b , yields 



= f\ h f eaa f/ c + f ebd f d ca + f ecd f d ab = -tr (T g T c T e ) - tr (T c T e T 9 ) + C 2 f ecg . (C.7) 
Using the cyclic invariance of the trace, this yields the "triangle rule" 

tr (T 9 T c T e ) = ^-C 2 f gce . (C.8) 

Eq. (C.8) is graphically depicted in (4.10). 

Invariant 4-index tensor d abcd : Next we consider the following totally symmetrized trace of four (adjoint) 
representation matrices 

d abcd ._ tr ^ T {a T b T c T d}^ ^ (C 9) 

which is (/-invariant by construction. This invariant arises when considering traces of four matrices T, Eq. (C.2). The 
result is given in (C.12) below. To derive it, observe that for traces of more than three generators T, which cannot be 
reduced using (C.8), one can permute two T's, with the aim of creating a loop of 3 with the remaining T's, which, in 
turn, can then be reduced using (C.8). For a trace of four T's, this reads 

tr (T*T b T c T d ^J = -\c 2 f ahh f cdh + tr (r b T a T c T d ^j . (CIO) 

This also tells us that 

^rpCLrpbrpCrpd^ ^ ^rpb rpCL rpd rpC^ ^ ^rpdrpCrpbrpO^ ^ 

where the second relation is obtained using the cyclic invariance of the trace. We now want to calculate a general trace 
of four T's. First, by using (C.ll), and the cyclic invariance of the trace, we find that of the 6 possible permutations, 
which leave the first index unchanged only 3 are independent. These are K x = tr (r a T b T c T d ) = tr (r*T d T c T b ), 
K 2 = tr (T a T c T d T b ) = tr (r a T b T d T c ), and K 3 = tr (T a T d T b T c ) = tr (r a T c T b T d ). The totally symmetrized 
trace, defined in (C.9), can now be expressed in terms of the Ki as: d abcd = \{K\ + K 2 + K%). Writing K\ = 
d abcd + | [(Ki — K 2 ) + (Ki — K 3 )], and using (C.6), we can rewrite each of these terms with the help of d abcd and 
f's as 

tr ^J ,ar J ,br J ,cr J ,c ^j jabcd _|_ ^2 j-adh jbch _ jabh jco 



(C.12) 



The invariant d 2 is now defined by 

J_ d abcd d abcd = 9l + d 2 (C.13) 
A/ a d 24 

where A/ad is the dimension of the adjoint representation. Note that (C.12), (C.13), (C.8) imply 

tr (r a T b T c T d ^j tr (r a T b T c T d ^j = N ad + d 2 ^j (C.14) 

The l.h.s. can graphically be viewed as the "cube"-invariant, discussed in section 2, and depicted in Fig. 2.3 of the 
same section [recall (C.2)]. 

In section C.2, we show that for SU(iV) 

d 2 = ^N 2 . (C.15) 

(The quadratic Casimir is C 2 = N in our conventions.) This is in agreement with the results of Ref. [22]. Hence, for 
SU(iV), the term d 2 in (C.13) is subleading in N, as compared to the first term. This subleading iV-dependence of d 2 
is also true for all the remaining classical groups, which follows from (C.13) and (C.27). 
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C.2 SU(iV) 

In this section we present a derivation of the value of the invariant c?2 for Q = SU(./V), i.e. of (C.15), which provides 
an independent check of this result given in Ref. [22] 15 . 

We start by recalling the generators in the (complexified) Lie algebra of SU(iV) in the fundamental representation 

X a = X? := {matrix with 1 in column a, row a; elsewhere} , (a, a = 1, N) (C.16) 

where the adjoint index a = {%}. These satisfy 



8%X£-6%X%, (C.17) 



which yields the structure constants in this basis 



f* = = - • (C18) 

The Killing form is given by 

ab _ a/3 ._ ^1 / T arpb\ — _( f, ac f b A — _9 [ Ffi$- - 

^ iv V )~ N V d Jc )~ \°r a N p ) 

= —2 ( projector onto the adjoint) (C.19) 



and its inverse 



One easily finds 



Vab = rj af3 = — iOZO f} -—o a o- /3 



VabV = (N — 1)(N + V) = 7V a d = dimension of adjoint representation . (C.21) 

Since we use r\ to raise and lower indices a = {%} of structure constants which are traceless [see (C.18)], one can 
also use the simplified form 

-1 



Vab v ab - -y * ? 



■ smlp = -ftifi (C- 22 ) 

instead of r) ab , for calculational convenience. Writing (T a ) c b = f c ab we obtain 16 : 

_l_ tt ^ = _ N (C 23) 

_L tt (^T 6 ) tr (T a 'T v ) r, aa , Vw = N 2 (C.24) 

_l_ tr tt ^ T a' T 6' TC '^ %a , wr?cc , = 1^3 (C 25) 



l_ tr ^ arbTCT ^ tj . ^^c'^ Vaa'riwVce'Vdd' = \n 4 + ^N 2 (C.26) 



Comparison of (C.26) with (C.14) yields d 2 = l^ 2 , in agreement with (C.14), and (C.15). 



15 This method can also be used to calculate higher invariants [36]. 
16 Again we use a computer to do the algebra. 
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C.3 Other groups 



Besides SU(iV) we also consider SO(N) and SP(N). The results of Ref. [37] yield: 



SU(iV) : 


c 2 


= N 


^2 


= 3^2 
2 JV 


SO(iV) : 


c 2 


= N- 2 




= 24 - § TV + 6/V 2 - §JV 3 


SP(iV) : 


c 2 


2 


<* 2 


= 3 , 45 M , 3 jy2 , _3_ iy3 
2 ^ 32 JV ^ 8 JV ^ 128 JV 



(C.27) 



We have already quoted these values for the group theoretical invariants on figure 2.1, but repeat them here for the 
convenience of the reader. Note that one can always normalize the 1-loop coefficient in the (3-f unction for g [the 
term oc g 2 in (2.3)] to 1/2, by rescaling g by a constant. This means that the normalization-invariant quantity which 
enters at 4-loop order is c^/Cf- This allows us to perform the following checks on (C.27), by using well known 
isomorphism between the corresponding Lie algebras: 
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3[32 + iV(14 + iV)] 
8(2 + TV) 3 



(C.28) 
(C.29) 
(C.30) 
(C.31) 



D Some elementary Integrals 

In this appendix, we consider some elementary integrals, quoted in the main text. Consider two points z a and in 
the complex plane, which are well inside a circle of (large) radius R centered at the origin. It is then elementary to 
establish the following result: 



d 2 3 



-7rln \zb — z a \ 2 + ir In R 2 + tt In 



\i\<R ~ Za)($* ~ Z%) 

Furthermore, for |j — z a \ < a < \z a — zj,\ 

h-z a \<a~ ° (3 " z a )(r ~ z* b ) 
Finally, this implies upon taking the limit of R — > oo, 
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/ d 2 3 

J\i~ z a\,\i-Zb\>a 



(3 - Z a )(l* ~ Zl) \l~Z a \ 2 

as long as \z a — Zb\ > 2a (up to terms of order a 2 which are neglected). 
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